本章开始我们介绍一下Smartseq与10X的数据结构
Smartseq首先,我们可以从网上下载示例数据,进行R里面的相关分析
wget https://ndownloader.figshare.com/files/10038307 unzip 10038307 wget https://ndownloader.figshare.com/files/10038310 mv 10038310 FACS_metadata.csv wget https://ndownloader.figshare.com/files/10039267 mv 10039267 FACS_annotations.csv对于Smartseq的数据,我们先读入FACS/Kidney-counts.csv这个文件
dat = read.delim("FACS/Kidney-counts.csv", sep=",", header=TRUE) #查看前5行5列dat[1:5,1:5]dat[1:5,1:5]如上图所示,第一列是基因的名字,而列名是细胞的谱系我们要将这个表达谱做成一个矩阵形式,那么我们把第一列基因名做成行名,并把第一列删除
rownames(dat)